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Abstract 

The Schrodinger equation with attractive 5 potential has been pre- 
viously studied in the supersymmetric quantum mechanical approach 
by a number of authors, but they all used only the particular superpo- 
tential solution. Here, we introduce a one-parameter family of strictly 
isospectral attractive <5 function potentials, which is based on the gen- 
eral superpotential (general Riccati) solution, we study the problem in 
some detail and suggest possible applications. 

The 5(x) (pseudo)potential is a well known 'zero-range' potential with 
applications in solid state physics (lj and many other areas. It has been 
used as a textbook example for many mathematical procedures in quantum 
mechanics. One such technique, Witten's supersymmetric scheme ||, has 
been employed for the attractive delta potential by several authors H [|, ||]. 
However, in all those studies there is a missing point, namely all the authors 
so far used only the particular Witten superpotential Wo, which is related 
to the ground state wavefunction in the well-known way uq = e~ -f w ° , and 
no mention is made of the general superpotential, i.e., the general Riccati 
solution for the 5 potential case. In this work we present the supersymmetric 
approach to the attractive delta potential problem based on the general 
superpotential. 

To help the reader to better understand our problem we start with its 
underlying mathematical scheme. Thus, we consider a Riccati equation (RE) 
of the type W = —W 2 + V^cc) for which we suppose to know a particular 
solution Wo. Let W\ = Wo + u be the second solution. By substituting 
W\ in RE one gets the Bernoulli equation u = —u 2 — 2Wqu, which by 
means of u = 1/v is turned into the first order linear differential equation 
v —2Wqv — 1 = 0. The latter one can be solved by employing the integration 
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factor /o = e~ 2 f Wo , leading to the solution v = / " 1 (C + f x /o)> where C 
is an arbitrary integration constant. Coming back to the general Riccati 
solution, one gets 



Wi = W + 



/o 



W + 



dx 



ln(C + | X /o) 



The point now is that in the process of factorizing the one-dimensional 
Schrodinger operator —d 2 /dx 2 + V\{x) the aforementioned Riccati solutions 
occur in the non-operatorial part of the factorization operators as follows. 
Wq occurs in the case of Witten's factorization (—d/dx + Wo)(d/dx + 
Wo){= AqAq), whereas W\ occurs for Mielnik's factorization || (—d/dx + 
Wi){d/dx + W x )(= A{Ai). Notice that [Al,A ] = 2W' Q , whereas [A\, Ax] = 
2W 1 . We further notice that vTo is the ground state (nodeless) wavefunc- 
tion of V\ and AVb r 

potential V\, leading to a new potential V\ do 



2Wq is the Darboux transform contribution to the 

Vi - 2Wn = V 2 , which in 



supersymmetric quantum mechanics is known as the supersymmetric part- 



ner of the initial potential V\. Even more interesting is that 



can 



be interpreted as the ground state wavefunction corresponding to Mielnik's 
superpotential (see below), and AVi = —2W 1 can be thought of as the 
general Darboux transform part in the potential. Therefore, there is a one- 
parameter family of Darboux potentials given by V\ t Di = V\ — 2Wi, which 
are strictly isospectral to the initial one, in the sense that each member of 
the family has the same supersymmetric partner V2 and the same energy 
eigenvalues and scattering amplitudes as V\. In terms of the ground state 
wavefunction of V\, ifto = yfo, each member of the strictly isospectral family 
of potentials reads 

V iso ,i = Vi + AVi = V 1 (x)-2-^ In (Ci + f /o) (2) 



or 



Vi(x) 



4^0 V'o 



+ 



2^ 4 



(3) 



For all half-line potentials the lower limit of the integral term is zero, whereas 
for the full-line potentials is —00. The ground state wavefunctions of this 
family are of the type ipo : i s „ = — . %°x . Indeed, one can write 



Wi 



— In 

dx 



— lnifjo,iso , 
dx 



(4) 
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which is the supersymmetric formula introducing the superpotential in terms 
of the ground state wavefunction. If one consider these isospectral functions 
as quantum mechanical wavefunctions, the problem of the normalization 
constant should be contemplated. It is easy to see that the normalization 
constant is N iso = \/C(C + 1) Q and as such C is not allowed to be in 
[—1,0]. The C = limit is known as the Pursey limit ||, whereas the 
C = — 1 limit is the Abraham- Moses limit M. However, in the present work 
we shall consider both the case with the normalization constant included 
and the case without it. 

Let us now pass to the attractive g6(x) potential, where g < gives 
the strength of the interaction (the binding power). It has been shown 
that Wq = |sign(x) In other words, Aq = d/dx + |sign(x) and A = 
— d/dx + |sign(x). Indeed, one cannot use the Heaviside step function as the 
superpotential since its square is not a constant. Therefore, one should work 
with the sign function, which is a combination of step functions. A^tpo = 
implies -00 = \/—g/2e 9 ^ x ^ 2 and this ground state wavefunction is the only 
one of the bound spectrum at the energy Eq = —g 2 /4. Thus, this state 
will be deleted from the spectrum of the partner potential, which is purely 
repulsive. However, the situation is by far the more interesting in the case 
of the strictly isospectral construction as one can see in the following. 

A simple calculation shows that 

l(x) = J* ^{x')dx> = -isign(x)effW + + \ . (5) 

Thus one gets 

t/ x r \ v \ , o 2 Csign(x)e-fl x l 

Visa = gSiso(x) = gS{x) + 2g — — _ , , (6) 

(1 — Csign(x)e 9\ x \y 



and the isospectral wavefunction reads 



sign(x)e" 



-91*1/2 



feo = - V^2~g ^C(C + 1) ^ , _ glxl , , (7) 
v (1 — Csign(x)e 9\ x \) 

where C = 2C + sign(x) + 1. The eigenvalue corresponding to the isospec- 
tral wavefunction is the same as for the common delta bound state, i.e., 
Eq = —g 2 /4. The analysis of Eqs. (6) and (7) shows that possible singu- 
larities are to be found for C in the interval (-1,-1/2], which is excluded 
when one considers normalizable isospectral wavefunctions. However for 
non-normalizable solutions these singularities should be taken into account. 
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The plots we did for the isospectral potentials as a function of the isospectral 
parameter (figure 1) display a shallow potential well on the negative half- line 
moving toward the origin where it is absorbed by the delta singularity there, 
and on the positive half-line a tail dying off at increasing C We also present 
plots showing the behaviour of the normalized isospectral wavefunctions for 
the same values of the C parameter as for the potentials (see figure 2). 
Moreover, figures 3 and 4 display the moving singularity structure when we 
do not introduce the normalization constant in equation (7). In summary, 
we believe that the strictly isospectral extension of the attractive 5 poten- 
tial introduced here may be relevant for many applications, once one allows 
for a physical origin of the C-dependence. For example, the parameter C 
may express the effect of static and/or moving distant boundaries, as well 
as sample-size dependence Jl^, [□]]. If one does not discard as unphysical 
the non-normalizable isospectral solutions, one may think of the isospectral 
method as allowing to introduce singularities in both wavefunctions and po- 
tentials which apparently are required to explain the extra losses of ultracold 
neutrons at the walls [12 1 . 
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Figure Captions 
Figure 1. 

Darboux potential contributions for C equal to .00001, .10001, 1.10001, 
5.10001, for g = —1 from left to right. 
Figure 2. 

The corresponding isospectral wavefunctions for the same values of C 
as in figure 1 showing how in the infinite limit of C one recovers the orig- 
inal 5 wavefunction. Actually, for rather low values of C, the isospectral 5 
wavefunction is already very close in shape to the original one. 

Figure 3. 

Darboux potential contributions for C equal to -1.4, -0.9, (up), and, -0.6, 
-0.3 (down) for g = — 1. 
Figure 4. 

Non-normalizable isospectral wavefunctions for the same values of C as 
in figure 3, together with the original ground state 5 wavefunction displayed 
in the first plot of the figure (g = —1). 
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